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Abstract 

We report the solution of the persistent current in two coupled rings in the presence of external 
magnetic fluxes. We showed that the magnetic fluxes modify the global phase of the electronic 
wave function for multiple connected geometry formed by the coupled rings. We obtained an exact 
solution for the persistent current and investigated the exact solution numerically. For two large 
coupled rings with equal fluxes, we found that the persistent current in the two coupled rings is in 
fact equal to that in a single ring. This theory explains the experimental results observed in a line 
of sixteen coupled rings. (Phys. Rev. Lett. 86, 3124 (2001).) 
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Introduction. In Quantum Mechanics, the wave function is path dependent and is 
sensitive to the presence of a vector potential caused by an external magnetic flux. In a 
closed geometrical structure, such as a ring, the wave function is changed by a measurable 
phase [1], causing all the physical properties to become periodic [2]. When a mesoscopic ring 
of normal metal is pierced by a magnetic flux $ [2], the boundary conditions are modified, 
leading to a famous theorem of periodic properties with the flux period $o = h/e and to a 
remarkable phenomena [3] of a non-dissipative persistent current [3-5]. 

One way to classify the closed geometrical structure is by using the number of holes 
formed on the closed geometrical structure. For a closed surface, the number of holes 
formed thereon is often referred to as a genus number g [6, 9]. For example, a genus number 
g = 1 describes an Aharonov-Bohm ring geometry, while a genus number g = 2 describes 
two rings perfectly glued at one point to form a character "8" structure. 

In spite of the predictions that the global geometry [1] has to affect the wave function, 
persistent currents were believed to be a specific property of isolated single rings. Only a 
recent experiment of 16 GaAs/GaAlAs connected rings [12] has demonstrated the possibility 
of persistent current in connected rings, namely, in a material with a large genus number. 
The experimental result showed that the current in a 16-rings structure is approximately 
the same as the current in a single ring, which disagrees with the existing theoretical results 
[11]. 

In this Letter, we report that for multiple connected geometries, such as a geometry 
with two or more holes (or connected rings), the geometry modifies the global properties 
of the wave function, and the presence of magnetic fluxes generates persistent currents 
with complicated periods. We present an exact analytical solution for the eigenvalues and 
compute the persistent current for two coupled rings with a character "8" structure for 
two different fluxes. We solve the problem by modeling the gluing of the two rings using 
fermionic constraints [14] with anti-commuting Lagrange multipliers, which can be viewed 
as a resonant impurity strongly coupled to the two rings [13]. 

The analytical results are investigated numerically. When the two fluxes on both rings are 
the same, we find a simple relation between the single ring (g = 1) current, I^ 9=1 \flux\ N), 
and the double ring (g = 2) current, I^ 9=2 \flux; N). At T = 0, we define I^ 9=2 \flux] N) = 
r(N)I^ 9=l \flux] N), where r(N) is a ratio between the two currents. The ratio r(N) is a 
function of the number of sites N and obeys r(N) — > 1 for N — > oo. This result explains 
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the experimental puzzle for a line of 16 GaAs/GaAlAs connected rings [12]. 

Exact Solution for Two Coupled Rings. We consider the Hamiltonian H for two 
spinless fermionic rings in the absence of a magnetic flux. The rings obey periodic boundary 
conditions. For each ring, the point x is identified with the point x + L. The two coupled 
rings with the character "8" structure (i.e. g = 2) are obtained by identifying the middle 
point x = L/2 of the first ring with point x = of the second ring, i.e. C\{L/2) = 6*2(0) 
and C 1 f (L/2) = C^(0). This identification is equivalent to two fermionic constraints, 
Q = Ci(L/2) — C 2 (0) and Q + = (L/2) — (7^(0). Since the constraints are fermionic, they 
can be enforced by using anti-commuting Lagrange multipliers, /i + and //. Following ref. 
[14], we introduce the Hamiltonian with the constraints, 

H T = H + fi+Q + Q+/x 

The unusual physical meaning of the anti-commuting Lagrange multipliers can be viewed 
as a fermionic impurity [13], which mediates the hopping of the electrons between the two 
rings. The two rings of length L are threaded by a magnetic flux $ a , where a — 1,2 
(for each ring). In order to observe the changes of the constraints in the presence of the 
external flux, we perform the following steps. In the absence of the external flux $ a , the 
annihilation and creation fermion operators obey periodic boundary conditions C a (x) = 
C a (x + L) and C+(x) = C+(x + L), where a = 1,2. The genus g = 2 is implemented 
by the Fermionic constraints Q = C^L/2) - C 2 (0) and Q+ = C?(L/2) - C 2 + (0), and 
Hamiltonian H = —tJ2 2 a=1 Y^x=G^ a \Pa{ x )^a{ x + °) + h.c.]. The length of each ring is 
L = N s a, where N s is the number of sites and a is the lattice spacing. When the external 
magnetic flux $ a is applied the Hamiltonian Hq is replaced by H. The Hamiltonian H is 
obtained by the transformation C a (x) —> exp[i-^ J* A(x')dx']C a (x) = ip a ( x ) an d C£(x) — > 
C£(x) exp[— ij^ Jq A(x')dx'\ = ^(x). Here A(x) is the tangential component of vector 
potential on the ring. The relation between the flux and the vector potential on each ring 
is j- Jq A(x)dx ~ j-A a N s a. 

The flux $ a on each ring a = 1, 2 gives rise to a change in the boundary conditions, 
ip a {x + N s a) = 4> a (x)e i<Pa and + A^ s a) = ^+(x)e _i *' a , where ^ a = 2tt(^) = 2tt|^ = 
2ir<p a . This boundary condition gives rise to a normal mode expression for each ring, ip a (x) = 
~7n ^u^q 1 e iK( ' n '' Pa ) ' x ip a (n) and a similar expression for ip^(x). The "momentum" K(n,ip a ) 
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is given by 

2tt 

K(n,<p a ) = -^(n + fi a ) 
N s a 

where n = 0, 1, ... ,N — 1 are integers with N = N s , and ip a = 2n0 a . In the momen- 
tum space, the Fermionic operators i/j a { n ) and ipt(m) obey anti-commutation relations, 
[i[> a (n), if)t(m)]+ = 8oi,p&n,m- The Hamiltonian for the two rings in the presence the external 
flux takes the form, 

(JV s -l)a N 3 —l 

H = -t^2 ^(^^(x + a) + h.c. = ^ ^ e(n,0 a )^(n)^ a (n) (1) 

a=l,2 x=0 a=l,2 n=0 

where e(n, (p a ) = — 2tcos[^(n + tp a )] are the eigenvalues for each ring. The Hamilto- 
nian in eq. [T] has to be solved together with the transformed constraints equations, Q = 
M^)e~ ilpl ~ MO) and Q+ = )e^ - ^ 2 + (0). 

The wave function for the genus g — 2 problem is given by the eigenstate \x) of the 
Hamiltonian in eq. HJ which in addition satisfies the equations Q\x) — and Q + \x) — 0. 
The constraint conditions are implemented with the help of the anti — commuting Lagrange 
multipliers /i and fi + . The Hamiltonian with the constraints takes the form, 

H T = H + fi + Q + QV (2) 

The Lagrange multiplier are determined by the condition that the constraints are satisfied 
at any time. Therefore, the time derivative satisfies the equation, Q\x) — Q + \x) = at 
any time. We will use the notations, [A, B] + = AB + BA and [A, B] = AB - BA. The 
Heisenberg equation of motion for the constraint Q is, 

ihQ = [Q, H T ] = [Q, H] + [Q, fi + Q + Q + fx] 

= [Q, H] + [Q, fi+] + Q - fi + [Q, Q} + + [Q, Q + } +f i - Q + [Q, //]+ 

= lQ,H} + [Q,Q+} + u (3) 

The rest of the anti-commutators in eq. [3] vanishes. The anti-commuting Lagrange mul- 
tipliers satisfy, [Q,fi + ] + = [Q,u] + = [Q + ,fi + ] + = [Q + ,u] + = 0. Since the constraints are 
fermionic, we obtain that they obey [Q,Q + ]+ = [Q + ,Q] + = 2. Therefore, the constraints 
are second class constraints [14]. From the condition Q\x) =0 and eq. [31 we determine the 
Lagrange multiplier field \x. 

fi=-[Q + ,Q]^[Q,H] = -^[Q,H] (4) 



fi + is obtained from the equation Q + \x) = 0, 

^ = [Q,Q + ]- + 1 [Q + ,H]= 1 -[Q + ,H] (5) 

The Hamiltonian Ht with the constraints and the Lagrange multipliers are used to 
compute the Heisenberg equation of motion for any Fermionic operator, O. (The Lagrange 
multipliers anti-commute with any Fermionic operator, i.e. [O, = [0,/i + ] + = 0.) 

ih^ = [6, H T ] = [6, H] + [6, //+Q] + [6, Q + ]u 

= [6, H] + [6, /i+]+Q - fi + [6, Q) + + [6, Q+]+n - Q + [6, fi] + 

= [d,H}-[d,Q} +f , + -[d,Q+}fM (6) 

We substitute in eq. [6] the solution for the Lagrange multiplier fields given by eqs. H] and [51 
We obtain a new equation of motion with a new commutator , which resemble the classical 
Dirac brackets [13]. 

ih^ = [0,H T ] = [6,H] - \6,Q + } + (\Q + ,Q} + )- l \Q,H} - [d,QU[Q,Q + ] + )- 1 [Q + , H] 
= [6,H] D (7) 

Eq. [7] shows that the Heisenberg equation of motion is governed by a new commutator, 
[O, H]jj. We will use this new commutator to compute the Heisenberg equations of motion 
for the creation and annihilation Fermionic operators ip a (x,t) and if)£(x,t), where a = 1,2. 

ihip a (x) = [ip a (x),H] D = [tp a (x),H] - ^[i! a (x),Q + ] + [Q,H] 

= -t[ip a (x + a) + ip a (x - a)] - -[8 a ,i8^ L / 2 e ltpi - 5 a , 2 4,o] 

i-t){e-^[M^ + a)+M^~ a)} + e**»[M^ + «) + U\ ~ «)]} (8) 

The ground state wave function is obtained from the one electron state, |% >= 
S Q =i 2 Si— o Zct{x)'ipa{ x )\® >> given in terms of the site amplitudes Z a (x). Using a 
normal mode momentum expansion, f a (n), i.e. Z a (x) = -j= J2n=o e lK ^ n '^ x f a {n), we find 
the following equations for the eigenvalues A and the amplitudes in the momentum space 
fa(n), 

(A - e(£ + <h))h{£) = ~ < n + ^) e ™/i(«) -^E e (« + ( 9 ) 

n=0 n=0 



and 



(A - e(£ + <h))f 2 (£) = ^J2< n + MM + ^J2< n + Vi)e im h{n). (10) 



2N ^ v r ,J y ' 2N 

n=0 ra=0 



We diagonalize eqs. I§1 and [TU1 by linear transformations, Sx((pi,X) = — ^{Lq 1 e (^ + 
0i)e %nt fi(£) and S , 2 (<^ , 2 ) A) = — ^fS^ e(£ + ^2) f2^) ■ As a result, we obtain the equation, 

( Si\ ( -{I 

M =0, where the matrix M is given by M = 

W 

we define Ai +) (<£ Q ,A) = Ai et,en) (<£ Q , A) + A% (<p a , A) and Ai '(^A) = Ai ei,en) (<£ Q , A) - 
Aa ((^ a , A), with the even and ogM representations given by, Aa ven ^ ((f a , A) = 

1 WiV-2)/2 £ (2m+y a ) 1 * (odd) / * ,v 1 WJV-2)/2 e (2m+l+g n ) w t 

2jvZ^m=0 A- £ (2m+0 Q ) allCl n « l^a, - 2JV L m =0 A- e (2m+l+^)' We COm P Ute 

det M = and obtain the characteristic polynomial which is used to determine the 
eigenvalues A. 

2[^r n) (0u A)A^ 2 , A) + A[ odd \^, X)^r n \02, A)] + [A< + >(fc, A) - A«(y3 2 , A)] = 1 

(11) 

Eq. [TT] is our main result for the genus g — 2 case. We observe that the matrix M is 
symmetric and the eigenvalues are real when the fluxes are equal, i.e. ifi = (f>2, or opposite, 
i.e. (p\ = —02- Fo r other cases, the eigenvalues can have imaginary parts, thereby giving 
rise to non-conducting states. 

Numerical Solution and Comparison with Experiment. We have numerically 
solved the secular equation [TU To compute the current, we sum over the current carried 
by each eigenvalue X(<fi,<f2) using the grand-canonical ensemble. The current in each ring 
a = 1, 2 is given by, 

i<r>(fc,fc)=- E ^.w ^^H ) (12) 

. ,* * . a^Pa Boltzman - 1 - 

where F( - w y V2> — / ^ m 'M = [1 + e K BoitzmanT l" 1 is the Fermi Dirac function which de- 
pends on the chemical potential Ef erm i and temperature T. The current is sensitive to the 
number of electrons being either even or odd. We use the grand-canonical ensemble and 
limit ourselves to a situation with even numbers of sites and a zero chemical potential, i.e. 
Efermi = (which corresponds to the half-filled case). In order to have a perfect particle- hole 
symmetry, we will restrict the analysis to the special series for the number of sites being 
N s = 2, 6, 10, 14, 18, ... , 2m + 2, where m — 0, 1, 2, 3 — For this case, we find that, when 




the fluxes are the same in both rings, the current for g = 2 has the same periodicity as that 
of a single ring, i.e. 

j(s=2)($ + $ ) = J(9= 2 )($). At temperatures T < 0.02 Kelvin, the line 
shape of the current as a function of the flux is of a sawtooth form (see figure 16). For other 
series N s ^ 2m + 2, the periodicity of the current is complicated. Using the experimental 
values given in the experiment [12], we estimate that the number of sites in our model should 
be in the range of N s = 50 ~ 150, the hopping constant should be t = ^ — hv fer™ — ~ q 
eV, and the temperature in the experiment should be T = 0.02 Kelvin. Using these units, 
we obtain that the persistent current is given in terms of a dimensionless current, I a (see 
figure 16 and figure lc) with the actual current value, Jq 9_2 ' ) = I a x 2.5 x 10~ 3 Ampere. 

a) Equal fluxes for g = 2: For this case the secular equation is simplified and takes the 
form of 4[A( even \0, X)A^ odd \p, A)] = 1. 

For N s = 2, we solve analytically the secular equation. We find that the eigenvalues are 
given by A(n, ip;N = 2) = r{N = 2)e(n,p>;N = 2), where e(n,<p,N = 2) = -2tcos[j^(n + 
<£>)], and n = 0, 1 are the single ring eigenvalues. The value for r(N = 2) is r(N = 2) = 
To find the eigenvalues for other number of sites, N s = 6, 10, 14, 18, 22, 26, 30, we numerically 
find the relation, A(n, ip; N) = r(N)e(n, tp; N), where n = 0, 1, . . . , N — 1 and e(n, p; N) = 
-2tcos[^{n + 0)} are the single ring eigenvalues. The function r(N) is given in figure la. 
This figure shows that the function r(N) reaches one for large N. Using the function r(N) 
given in figure la, we compute the current for the g = 2 case as a function of temperature. 

I^\ r ,N;T) = -'"f 1 ± H N).^N)]F C^-<;^^ z£te=J ) (13) 

^Tq aL P A Boltzman 1 

Figure lb represents the current computed from eqs. [TT] and [12] for N s = 30 sites 
at two temperatures T = 0.02 and T = 20 Kelvin. In this figure, the current is given 
in dimensionless units / plotted as a function of the dimensionless flux (p a = [—0.5, 0.5] 
(ip a = 27[0 a = [— 7T, 7r]). The solid line represents the single ring current and the dashed line 
represents the current for the genus g — 2 case. In figure 16, the ratio of the currents at 
T = 0.02 Kelvin is r(N = 30, T = 0.02) = 0.979. 

The experiment [12] was performed with 16 rings, using the condition that each ring 
has 50 sites, which gives an estimation of r = I w _ rings / I sing i e _ ring w [J (9=2) // S m 9 «e-rm(,] 4 = 
[r(T = 0.02, N s = 50)] 4 = [0.987] 4 = 0.95. The ratio of r = 0.95, is in the range of the 
experimental observation reported in ref.[12]. 

Next, we turn our attention to the values of the currents. The current reported in ref. 



12 is in the range of 0.5 nA. Our estimate for 30 sites given in figure 16 at T = 0.02 is 
/ = 10~ 3 , which corresponds to a current of I^= 2 ) = / x 2.5 x 10~ 3 v4 = 2.5 x 10 3 nA. 
This means that our current is 10 3 larger then the reported current [12]. The origin of this 
discrepancy might be due to the following factors: the effective number of sites is larger then 
30; the temperature in the experiment might be higher then T = 0.02 Kelvin] and the elastic 
scattering length leiastic in the experiment is smaller than the length of the ring L. Assuming 
a strong 2Kp impurity scattering in the presence of the electron-electron interaction can 
cause a significant current suppression [7,8]. For an impurity with a transmission coefficient 
t < 1 and a repulsive electron-electron interaction with a Luttinger parameter K c < 1 the 
reduction of the current is (|£|2)(^ini)2(^ c - 1 -i) ) [ 7)8 ]_ The ratio (^«) ~ _i_ wit h K c = 0.6 

gives a suppression of 1000 for the current, in agreement with the observed currents. 

At T = 20 Kelvin, the value of the currents are in the range of 7 nA and the reduction 
of the current is larger in comparison with the T = 0.02 Kevin. 

We have checked the behavior of the current as a function of the number of sites (for 
the single ring and the double ring) at a fixed flux 0.001 and different temperatures. At 
T = 0.02 Kelvin the current scales like (for the double ring and the single ring) up to 
Ns ~ 60 . Increasing the temperature to T = 1. Kelvin the current decreases faster than 
-j- for Ns > 30 . This results confirm that at low temperature the Persistent current in a 
double ring decreases linearly with the length. 

b) Two coupled rings with opposite fluxes, i.e. (p — <f% — —(f>2- 

For N s = 2, the eigenvalues are the same as the one obtained for the same flux 
case. For N s = 6, 10, 14, ... , 2m + 2, we solve the secular equation given in eq. dU 
and compute the eigenvalues. In figure 26, we plot the total energy as a function of 
the opposite fluxes at T = 0.027^ for 30 sites, E^ =2 \-0,0,N s = 30, T = 0.02, K) = 
V™~ ( ^~ 1 [A(— <^>, 0)F( ( ' X( ' jf'^ E f £ ™^ \\ . We observe that the total energy has a chaotic struc- 
ture due to the strong Backscattering caused by the common point of the two rings. For 
comparison, we also show in figure 2a the total energy for equal fluxes. The energy is 
parabolic for small fluxes and the current is proportional to the flux. For large values of flux 
the total energy E^ 9=2 ^ is a periodic function with the period $o- 

c) We have also considered the case with unequal fluxes. In ring one, the flux was fixed at 
the values (p\ = 0.1, 0.2, 0.3, 0.4, and in ring two, the flux <p 2 changes continuous from —0.05 
to 0.05. For this case, the current in ring two is linear in the flux 02 and is independent on 
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the flux in ring one, (p\. 

d) An effective four terminal circuit model, which explains the physics of the g = 2 
coupled rings. 

The persistent current problem is mapped into a closed circuit with a voltage source V. 
We use the relation between the flux and the voltage [7], ^ = j^f/ermi( From this 
relation, we find that a perfect ring is equivalent to one conducting channel connected to a 

2 

battery with a voltage V. This gives a current I — ^V. The g — 2 case is described by 
two circuits coupled at x = L/2. The first ring with a flux is equivalent to a circuit that 
start at x = runs to the common point x = L/2 and back to x — L. We attach a voltage 
source V\ = (| ) _1 ;^f/ e rmi(fr;) between the points a; = and x = L. We do the same for 
the second ring, where we attach a voltage V2 = —{^)~ l j^ v fermi{i^)- 

At the coupling region, we split the point x = L/2 into two points O and O', such that a 
/wg/i resistance is attached between the points O and O'. Now, we consider separately 
the equal and opposite fluxes. For equal fluxes, we have to use two opposite voltages 
Vi = — V 2 = V/2. For this case, the current does not pass trough the high resistance region 
between the points O and O'. Instead, it forms a loop with twice the length of the single 
ring with a total doubled flux and an effective voltage V. Therefore, the current is the same 

2 

as in a single ring, i.e. I = \V ■ 

Next, we consider the situation for opposite fluxes. In this case, the two batteries obey 
V± — V2 — V/2. As a result, the current has to flow through the high resistance region O 
to O'. The currents are opposite in each ring and their value is determined by the high 
resistance i?^. We find, / = — 2 — 0. 

Summary. In this Letter, we have introduced a method which solves the problem of 
the global phase of the wave function for geometrical structures with holes, i.e. high genus 
materials. This method is applicable to a variety of mesoscopic systems where coherency of 
wave function is important. 

We have found an exact solution for the persistent current in two coupled rings. By 
numerical calculations, we have computed the current dependence on the flux, temperature, 
and the number of sites. 

This theory resolved the experimental puzzle [12] that the persistent current in many 
coupled rings is the same as that of the single ring. 
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FIG. 1: (a) The ratio of the double to single ring currents I(g = 2;N)/I(g = 1;N) = r(N); (b) 
The single ring (solid line) and the double ring (dashed line) currents for N s = 30 at T = 0.02 
Kelvin; and (c) The single ring (solid line) and the double ring (dashed line) currents for N s = 30 
at T = 20.0 Kelvin. 
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FIG. 2: (a) The total energy for equal fluxes, / = <pi = <p 2 E^= 2 \f,f;N s = 30, T = 0.02K); 
and (b) The total energy for opposite fluxes, / = — = </?2 for 30 sites at T = 0.02 Kelvin 
E(9=2) (-/,/; iV a = 30,T = 0.02K) . 
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